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An automatic pipeline analysing solar-like oscillating targets tested on 
CoRoT and simulated data 
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The launch of the Kepler mission on 7th March 2009 opened a new bright future for the search of extra-solar planets while 
a huge amount of stars will be observed leading to the opportunity to better understand stellar evolution. This will allow us 
to probe different regions in the HR diagram and put more constraints on the stellar models. Up to now the asteroseismic 
missions such as MOST and CoRoT were providing some solar-like stars at a very slow cadence. But to study the several 
hundreds of solar-like oscillating stars that will be observed during the Kepler survey phase, an analysis devoted to one 
star is impossible if we want to have as much information as we can in a small period of time. We describe here our 
pipeline, A2Z, which calculates the global parameters of the stars (rotation period, mean large spacing of the acoustic 
modes, maximum amplitude of the modes), fits the modes globally, and estimates the radius and mass of the stars. This 
pipeline has been tested on simulated stars and applied to real data from CoRoT. 
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1 Introduction 

The study of stellar evolution is going through an important 
revolution. Helioseismology had already allowed us to im- 
prove our knowledge on the solar interior, such as density, 
rotation, position of the base of the convection zone, and He 
content, showing the importance of this useful tool. 

During the last decade, the development of asteroseis- 
mology started in many ways: through ground-based ob- 
servations campaigns and thanks to many space missions 
allowing the observations of solar-like oscillations over a 
longer period of time (WIRE, MOST, CoRoT). 

More recently, the launch of the NASA mission, Kepler, 
dedicated to the search for exoplanets has also extended the 
study of stellar oscillations to a few hundreds of solar-like 
oscillating stars. 

To face the arrival of a huge amount of data, several 
teams have developed different pipelines (Huber et al, 2009; 
Hekker et al. 2010) in the framework of the AsteroFLAG 
group. We will describe the A2Z pipeline (Mathur et al. 
2010a). This pipeline is divided into several packages: the 
search for the frequency -range of the p-mode bump by look- 
ing for the large spacing (average distance in frequency be- 
tween consecutive nodes for the same mode degree), the 
fitting of the background, the estimation of the maximum 
amplitude per radial mode as well as the central maximum 
frequency. At this point, we can try to infer the radius and 
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the mass of the stars from all this information (e. g. Stello et 
al. 2009). Then, if the signal-to-noise ratio is high enough 
we can obtain the characteristics of the p modes through a 
global fitting on the power spectrum. For the first package, 
we also have several methods that enables us to cross-check 
the results and to attribute a confidence level. After describ- 
ing the methodology used in the pipeline, we will apply it 
to the Sun and to a few CoRoT targets. 

2 Pipeline packages 

The A2Z pipeline is constituted of 3 parts: the first one for 
retrieving the global parameters of the acoustic modes as 
well as a first estimation of the mass and radius using the 
scaling laws, the second one to fit the modes globally, and 
the last one to estimate the radius and mass of the stars by 
using stellar modelling. We will focus here on the determi- 
nation of the global parameters of the acoustic modes, the 
estimation of the guesses, and the inference of the radius of 
a star. 

If the signal-to-noise ratio is high enough the global- 
fitting code can be run and it will infer the central frequency 
(ttj), the amplitude (A{), the width (r^), and the rotational 
splittings (8ui) of each mode. 

To compute the Power Spectrum Density (PSD), we use 
the Lomb-Scargle algorithm, which takes into account the 
irregularity of the sampling rate of the time-series. 

In the case that there are gaps in the time series, we use 
the in-painting algorithm described by Sato et al. (2010). 
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Fig. 1 Workflow diagram of the A2Z pipeline. The circles are the inputs and outputs of the different packages (boxes). 
Each output parameter is given with its uncertainty (extracted from Mathur et al. 2010a). 



For very low signal-to noise ratio data, we also plan to apply 
the curvelets (Garcia et al. 2010). 

2.1 Average Photospheric Rotation 

Some information on the surface rotation of the star can give 
some constraints on the rotational splitting and the inclina- 
tion while fitting the modes. Indeed, it is very difficult to 
estimate these two values only with the fitting procedure as 
they are mostly correlated because of the short lifetime of 
the modes (Ballot et al. 2006, 2008) 

We apply two different methods to retrieve the rotation 
period of the stellar surface. First we can compute the PSD 
and look at the highest peak at low frequency, below 5 /iHz. 
However, sometimes, this method does not enable us to dis- 
tinguish between the fundamental period and its first har- 
monic. Thus we also use a time-frequency analysis, with 
the wavelet tools. We use the Morlet wavelet to calculate 
the correlation between this wavelet that we slide along the 
time-series. The frequency of the wavelet takes many values 
in a given range. Fig. [2] shows the results, wavelet power 
spectrum (WPS) for 1 month of solar data taken by the 
VIRGO (Frohlich et al. 1995) instrument aboard SoHO. If 
we collapse the WPS along time, we obtain the global wave- 
let power spectrum (GWPS) on the right panel of Fig. [2] We 
see on Fig. [2] a large dark pattern showing that there is some 
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Fig. 2 Wavelet power spectrum for 1 month of VIRGO 
data (left). Global power spectrum (right) (from Mathur et 
al. 2010). 

continuous power along time around a period of 27 days 
as well as a large peak in the GWPS. We can also provide a 
confidence level (95% is the value chose in our study). Thus 
we can give an estimation of the rotation period of the solar 
surface with more than 95% confidence level. 

2.2 P-Mode Range Search 

We are looking for the frequency range of the p-mode ex- 
cess power as well as an estimation of the mean large spac- 
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ing, (Av). We use two different methods to cross-check our 
results and to select only those for which we are confident. 
Method 1: we estimate (A^) by calculating the Power Spec- 
trum (PS) of PS (from 200 to 6000 Hz). We compute the 
PS of PS (or PS2) on sub-series of 600 /iHz. We obtain 
the frequency-range of the p modes, f m %n an d fmax with 
a given confidence level: from 85 to 95% (level 2), from 70 
to 85% (level 1), and below 70% (level 0) (see Fig. ID. 
Method 2: PSSPS (Power Spectrum of a Short or Small part 
of the Power Spectrum). Here, we estimate (A^) from the 
standard stellar parameters. Then, we search for the spacing 
in the PS of a short slice of the PS where the p-mode excess 
power appears (Regulo & Roca Cortes 2002). 




2000 4000 6000 6000 

Frequency (^wHz) 



Fig. 3 Frequency -range of the p modes for 1 month of 
VIRGO data. 

We have a 3-step strategy, which is the following: 

- Run method 1 and 2; 

- Run method 2 with the results of method 1 as inputs; 

- Run method 1 with the results of method 2 as inputs. 

To select the most reliable results, we dismiss the step in 
which we obtain the highest number of common stars with 
level of 0. Among the two remaining steps in which we re- 
move the stars with the confidence level of 0, we choose the 
step where we have the highest number of coincidences. 

2.3 Background Fitting 

For the background fit, we use a fit of three components, 
one or several Harvey laws (Harvey 1986), one power law 
for the activity effects, and a flat photon noise, W: 

B(y) = W+ * , a +av- h , (1) 
1 + (2v/v c ) 

where we have six free parameters. The PSD is fitted with- 
out the mode region and with a maximum likelihood esti- 
mator. 

2.4 Bolometric Amplitude of the £=0 mode 

We look for the amplitude envelope of the PS. First, we 
smooth the PSD corrected from background over 1 or 2 



x(A^). We fit the PSD with a Gaussian to get the maxi- 
mum power Pmax and v max . We convert the power into the 
amplitude per radial mode A max (Kjeldsen et al. 2008) and 
to bolometric amplitude (Michel et al. 2009). 

2.5 First Determination of the Radius and the Mass 

Based on the scaling laws (Kjeldsen & Bedding 1995), we 
can calculate the radius and the mass of the star with the 
following formulae: 

Rl _ ( 135 w*ww Tag y/ 2 

R Q \{Av)) V 3050 A 5777/ ( } 

M' _ / 135 T cB \3/2 

M Q \{Au)J V3050/ V5777/ 

2.6 Guessing Parameters 

We determine of the table of guesses of individual mode 
frequencies and variation of (Ai/) with frequency based on 
the results obtained with PMRS. We look for the highest 
peaks in the range [/ m i n , /max] and that are approximately 
separated by {Av), allowing an uncertainty of 15%. 

2.7 Radius Extractor 

Here, we use a pre-calculated grid of evolutionary models 
to obtain an initial guess of the parameters such as mass 
and age. With a minimization algorithm, we obtain a better 
estimation of the radius and mass. 

3 AsteroFLAG/Aarhus simulated stars 
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Fig. 4 Distribution of the correct results obtained with 
the simulated stars from the Aarhus-asteroFLAG simulator. 
Full black circles represent the stars for which the correct 
value of the large spacing was obtained. The squares are 
the stars selected to determine the radius: the converging 
models (green) and the diverging models (red). 

The A2Z pipeline was tested on 176 stars of magnitude from 
7 to 1 1 that had been simulated by the Aarhus-AsteroFLAG 
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Fit Radius (R 8 ) 

Fig. 5 The true radius versus the fitted input radius for a 
selection of the simulations. The continuous line is the x=y 
line. 
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Fig. 6 Smoothed bolometric amplitude per radial mode 
(rms) for the Sun, HD49933 (run 1 and 2), HD181906, 
HD181420, and HD49385. 



Table 1 Results of the A2Z pipeline on real data 



Star 


Period 


(Au) 


£ 


R 


e 




(days) 


(jzHz) 


OuHz) 


(R@) 


(Re) 


Sun 


27.8 


135.46 


3.26 


1.016 


0.023 


HD49933 


3.3 


86.14 


1.58 


1.398 


0.04 


HD181906 


2.8 


86.53 


1.69 


1.613 


0.032 


HD181420 


2.6 


75.35 


1.53 


1.42 


0.018 



simulator (Chaplin et al. 2008). These time-series are one- 
month long to simulate the data obtained with the Kepler 
mission. Following the 3-step strategy, we obtain a common 
mean large spacing between the 2 methods of the package 
PMRS for 71 stars, among which 67 are the correct ones. 
We have applied the other packages BF and BAL0 to cal- 
culate the background fit and the maximum amplitude per 
radial mode for these 67 stars. The amplitudes are often un- 
derestimated. Finally, the RadEx package was used on a few 
tens of stars to estimate their radii. Fig.|5]shows the compar- 
ison between the radius input and the result of the fit for a 
few tens of stars. 

4 CoRoT solar-like stars 

The A2Z pipeline has been applied to solar data from VIR- 
GO, to several CoRoT targets such as, HD49933 (Appour- 
chaux et al, 2008; Benomaret al. 2009), HD181906 (Garcia 
et al. 2009), HD181420 (Barban et al. 2009), HD175726 
(Mosser et al. 2009), HD49385 (Deheuvels et al. 2010), 
and HD 170987 (Mathur et al. 2010b), and to the Kepler 
data, solar-like oscillating stars and red giants (Bedding et 
al 2010; Chaplin et al, 2010; Stello et al. 2010). We give a 
few reseults in table 1. Fig. [6] shows the bolometric ampli- 
tude envelope per radial mode for the Sun and the CoRoT 
targets. 

Acknowledgements. This work has been partially supported by: 
the CNES/GOLF grant at the Service d'Astrophysique (CEA/Sa- 
clay) and the grant PNAyA2007-62650 from the Spanish National 
Research Plan. This work benefited from the support of the Inter- 
national Space Science Institute (ISSI), through a workshop pro- 



gramme award. It was also partly supported by the European He- 
lio- and Asteroseismology Network (HELAS), a major interna- 
tional collaboration funded by the European Commission's Sixth 
Framework Programme. SoHO is a space mission of international 
cooperation between ESA and NASA. CoRoT (Convection, Ro- 
tation, and planetary Transits) is a mini-satellite developed by the 
French Space agency CNES in collaboration with the Science Pro- 
grams of ESA, Austria, Belgium, Brazil, Germany, and Spain (Mi- 
chel et al. 2008). 



References 

Appourchaux, T, Michel, E., Auvergne, M., et al.: 2008, 
A&A 488, 705 

Ballot, I., Garcia, R.A., & Lambert, P.: 2006, A&A 369, 1281 
Ballot, I., Appourchaux, T, Toutain, T, et al.: 2008, A&A 486, 
867 

Barban, C, Deheuvels, S., Baudin, F, et al.: 2009, A&A 506, 51 
Bedding, T.R., Huber, D., Stello, D., et al.: 2010, ApJ, in press 

(arXivlOOl.0229) 
Benomar, O., Baudin, F, Campante, T.L., et al.: 2009, A&A 507, 

L13 

Chaplin, W.I., Appourchaux, T, Arentoft, T, et al.: 2008, AN 329, 
549 

Chaplin, W.J., Appourchaux, T, Elsworth, Y., et al.: 2010, Apl, in 

press darXiv: 1001 .05061 
Deheuvels, S., Bruntt, H., Michel, E., et al.: 2010, A&A, in press 

JarXiv: 1003 .4368t 
Frohlich, C, Romero, J., Roth, H., et al.: 1995, SolPh 162, 101 
Garcia, R.A., Regulo, C, Samadi, R., et al.: 2009, A&A 506, 41 
Garcia, R.A., Ballot, I., Mathur, S., et al: 2010, AN, this volume 
Hekker, S., Broomhall, A.-M., Chaplin, W.J., et al.: 2010, MN- 

RAS 402, 2049 

Huber, D., Stello, D., Bedding, T.R., et al.: 2009, CoAst 160, 74 
Kjeldsen, H., Bedding, T.R., Arentoft, T, et al.: 2008, Apl 682, 
1370 

Kjeldsen, H., and Bedding, T.R.: 1995, A&A 293, 87 

Mathur, S., Garcia, R.A., Regulo, C, et al.: 2010a, A&A 511, A46 

Mathur, S., Garcia, R.A., Catala, C, et al.: 2010b, A&A accepted 

(ai'Xiv 1004.4891) 
Michel, E., Baglin, A., Auvergne, M., et al.: 2008, Science 322, 

558 

Michel, E., Samadi, R., Baudin, F, et al.: 2009, A&A 495, 979 



© 2010 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinhcim 



www.an-journal.org 



Astron. Nachr. / AN (2010) 



793 



Mosser, B., Michel, E., Appourchaux, T., et al.: 2009, A&A 506, 
33 

Regulo, C, Roca Cortes, T.: 2002, A&A 396, 745 
Sato, K.H., Garcia, R.A., Pires, S., et al.: 2010, AN, this volume 
Stello, D., Chaplin, W.J., Bruntt, H., et al.: 2009, ApJ 700, 1589 
Stello, D., Basu, S., Bruntt, H., et al.: 2010, ApJ, in press 
(arXivlOOl.0026) 



www.an-journal.org 



© 2010 WILEY- VCH Verlag GmbH &Co. KGaA, Weinheim 



